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Abstract. Radial abundance gradients are a common feature of spiral galaxies, and in the case of the Galaxy 
both the magnitude of the gradients and their variations are among the most important constraints of chemical 
evolution models. Planetary nebulae (PN) are particularly interesting objects to study the gradients and their 
variations. Owing to their bright emission spectra, they can be observed even at large galactocentric distances, 
and the derived abundances are relatively accurate, with uncertainties of about 0.1 to 0.2 dex, particularly for the 
elements that are not synthesized in their progenitor stars. On the other hand, as the offspring of intermediate 
mass stars, with main sequence masses in the interval of 1 to 8 solar masses, they are representative of objects 
with a reasonable age span. In this paper, we present an estimate of the time variation of the O /H radial gradient 
in a sample containing over 200 nebulae with accurate abundances. Our results are consistent with a flattening of 
the O/H gradient roughly from —0.11 dex/kpc to —0.06 dex/kpc during the last 9 Gyr, or from —0.08 dex/kpc 
to —0.06 dex/kpc during the last 5 Gyr. 
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1. Introduction 

Radial abundance gradients in the Galaxy have been 
determined by a variety of objects, including HII re- 
gions, planetary nebulae (P N), B stars, o pen cl usters and 
cephe ids (H enry & Worthey |1999| , Maciel |2000| , Rolleston 



et al. 2000| ). Several elements have been investigated, es- 



pecially oxygen, sulphur, neon and argon in photoionized 



nebulae (Maciel & Quircza 1999 , Deharveng et al. 3000 ) 
and iron in open cluster stars and cepheids (Friel 1999 , 
Andrievsky et al. |2002a| , |2002bf ). In the case of the O/H 
ratio, recent determinations from HII regions, PN and 
B stars point to an average gradient in the range —0.04 
to —0.07 dex/kpc, taking into account an average uncer- 
tainty of 0.02 dex/kpc for most determinations. 

Abundance gradients play a distinctive role as a con- 
straint to chemical evolution models. In fact, several re- 
cently computed models, based on widely differing as- 
sumptions, point to the radial gradients as one of the most 
important constraints, especially when one takes into ac- 
count not only their magnitudes, but also their space and 
time variat ions (see for exam ple Ho u et al. 2000 . Chiappini 
and Alibes et al 



et al. 2001 



2001|) 



Planetary nebulae are particularly useful in the study 
of the gradients and their variations (Peimbert 1990, 



Maciel fc K oppen |1994[ Maciel |l997j , |2000j Peimbert & 
Carig i |199S| , Maciel & Quireza y])9|, Maciel & Costa 



[2002 ). They usually have bright emission spectra, so that 



they can be observed even at large heliocentric distances. 
Their derived abundances are relatively accurate, with un- 
certainties of about 0.1 to 0.2 dex, especially regarding 
the elements that are not synthesized by their progeni- 
tor stars, such as neon, argon and to some extent oxygen. 
On the other hand, as the offspring of intermediate mass 
stars, with main sequence masses roughly in the interval 
1 < M/Mq < 8, they are representative of objects span- 
ning a relatively large age interval, so that it is expected 
that groups of PN of different ages may display different 
abundance variations, thus reflecting the time evolution 
of the gradients. The PN distances are a possible source 
of error, as they are not as well determined as in the 
case of B stars or HII regions, for example. However, it 
has been shown that the use of different distance scales, 
both individual and statistical, apparently compensates 
for this uncertainty (see for example Maciel & Koppen 
1994). Moreover, it should be recalled that in order to 
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determine the PN gradients larger samples are generally 
used as compared with HII regions. 

In this paper, we present an estimate of the time vari- 
ation of the O/H radial abundance gradient based on a 
large sample of galactic PN with relatively accurate abun- 
dances and distances. From the observed oxygen abun- 
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dances, we determine the [Fe/H] metallicity using a cor- 
relation based on disk stars, and the progenitor ages are 
obtained through an age-metallicity relation. In section 2 
we describe our sample, in section 3 we discuss the [O/H] 
x [Fe/H] relation for disk stars. The ages of the objects 
are estimated in section 4, and we present our method to 
determine the abundance gradients as a function of age. 
Finally, Section 5 presents our results and discussion. 

2. The data 

The sample includes about 240 nebulae, most of which 
have oxygen abundances by number log (O/H ) and dis- 
tances from the sa mples of Maciel & Quireza ( 1999 ) and 
Maciel & Koppen ( 1994 ), to which the reader is referred 
for details on the abundances and references. About 40 
new nebulae have been included, basically from recent 
observations from our own group, secured at the 1.52m 
ESO telescope at La Silla and the 1.60 m LNA telescope 
in Brazil (Costa et al. |l996| , |l997t |2002[ ). Most of the new 
objects belong to a project to derive accurate and homo- 
geneous abundances of PN located near the anticentre di- 
rection. A detailed analysis of the new observational data, 
plasma diagnostics and abundances, as well as a study 
of the space variations of the abundance gradients along 
the galactic disk is given elsewhere (Costa et al. 2002| ). 
It should be mentioned that most O/H abundance deter- 
minations are based on the so-called "empirical" method, 
according to which the total abundances are obtained from 
ionic abundances with the use of ionization correction fac- 
tors for the unseen species. The ionic abundances depend 
on the plasma diagnostic parameters, namely the electron 
temperature and density, which are derived from selected 



1998, Torres- Peimbert & Peimbert 1997). Furthermore, 



line ratios (see for example Peimbert 1990 and references 
therein). This procedure implies average uncertainties of 
about 0.1 to 0.2 dex for the abundances of most elements 
and up to 0.02 dex/kpc for the derived gradients. These 
uncertainties are able to accomodate any systematic vari- 
ations with the ga lactoc entric distance (see for example 
Martins & Viegas 200C). In fact, the presence of abun- 
dance gradients can also be detected by considering abun- 
dances derived by other methods, such as photoionization 
models or detailed individual analysis, and the main dis- 
cussions refer to the magnitude of the gradients. 

Our sample consists only of disk nebulae, so that PN o f 
Types IV (halo objects) according to the Peimbert ( 1978 ) 
classifi cation system, or of Type V (bulge objects, see 
Maciel 1989 ) are not included. Most of the objects are 
classified as Type II (disk objects, intermediate mass pro- 
genitors), but PN of Type I (disk objects, large mass pro- 
genitors) and III (thick disk objects, showing large devia- 
tions from the galactic rotation curve) are also included, 
so that the average dispersion of the oxygen abundances 
of the whole sample is expected to be larger than in Maciel 



& Quireza (|199J). 

There are some suggestions that Type I PN may show 
some effects of ON cycling in their oxygen abundances, but 
current evidences are not conclusive (Peimbert & Carigi 



the amount of oxygen depletion from this process is ex- 
pected to be much lower than the average abundance 
dispersion and restricted to the objects of larger masses, 
which make up a small fraction of the sample. 

3. The [O/H] x [Fe/H] relation 

The determination of the radial O /H abundance gradient 
was p erformed in a similar way as in Maciel & Koppen 
( 1994| ), and simple linear fits have been obtained. The 
O /H abundances by number of atoms have been converted 
into the usual [O/H] abundances relative to the Sun us- 
ing the solar oxygen abundance log(0/H) R + 12 = 8.83 



(Grevesse & Sauval 1998| ), so that we have the relation 



[O/H] = log (O/H) - log (O/H), 



(1) 



Iron abundances are very difficult to determine from 
PN and their central stars, basically because the iron lines 
are very weak in PN spectra, and a presumably large 
fraction of iron is condensed in grains (see for example 
Perinotto et al. |1999 and Deetjen et al. 1999| ). In order 
to convert from [O/H] to the iron metallicities [Fe/H], we 
have used the extensive sample of disk stars by Edvardsson 
et al. ( 1993| ), from which we derive the relationship 



[Fe/H] = 0.0317+ 1.4168 [O/H] 



(2) 



valid approximately in the range — 1.0 < [Fe/H] < 0.5. The 
average uncertainty of the [Fe/H] metallicity is roughly 
0.1 dex. For this relation we have determined an average 
uncertainty for the slope of 0.0488, and for the intercept 
the uncertainty is 0.0158. The correlation coefficient is r = 
0.9536, which reflects the good fit of Eq. |2|, as can be seen 
in Fig. 0. 

In order to test the applicability of Eq. || to the galac- 
tic disk as well as to the solar neighbourhood, we have 
made some comparison with the [Fe/H] x [O/H] relation 
as obtained for different sets of galactic objects by dif- 
ferent authors. A conv enient set has been recently anal- 
ysed by Maciel ( 2002 ), in a study of the use of radial 
abundance gradients mainly derived from young objects 
in the determination of the [O/Fe] x [Fc/H] relation in 
the galactic disk. The main results are shown in Fig. |[ 
In this figure, the solid line shows the relation given by 
Eq. ||. The dotted curve s hows results of theoretical mod- 
els by Matteucci et a l. ( 1999| ), which follow models by 
Chiappini et al. ( 1997 ), and are representative of models 
predicting a [O/Fe] plateau for metallicities under solar, 
without a significant increase in the [O/Fe] ratio above 
0.5 dex for [Fe/H] < —2. The dashed line is the predicted 
relationship on the basis of the radial gradients of young 
objects in the galactic disk, at galactocentric distances 
4 < i?(kpc) < 12, as discussed by Maciel fl2002| ). The fig- 
ure also includes some representative observational data 



by Barbuy & Erdelyi-Mendes (1989), asterisks; Boesgaard 
et al. ( |l999l ), filled squares; Edvardsson et al. ([1993]), 
crosses; Israelian et al. ( 1998j ), solid dots; Spiesman & 
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Edvardsson et at. (1993) 




Fig. 1. Correlation between [Fe/H] and [O/H] obtained from the data of Edvardsson et al. (1993). 




[0/H] 



Fig. 2. A comparison of the average relation given by Eq. || (solid line), theoretical models by Matteucci et al. (1999, 
dotted line), the predicted relation derived from the gradients of young objects (Maciel 2002, dashed line) and selected 
observational data from a number of sources (see text). 



Wallerstein ( 1991 ), o pen ci rcles; Spite & Spite ( 1991 ), plus 
sig ns; Ta keda et al. ( 2002 ), empty squ ares; Mishenina et 



al. ( 2000 ), open stars, and Cavallo et al. (1997 ), filled stars, 
as recomputed by Mishenina et al. ( 2000 ). 

It can be seen from Fig. || that the observational data 
and models show a reasonable agreement, especially for 
metallicities close to and slightly lower than the solar 



value. For the lowest metallicities the spread is somewhat 
larger and the points lie, on the average, under the line cor- 
responding to Eq. ^, but the relatively sparse data are still 
reasonably well represented by the solid line. Part of the 
scatter in Fig. || may be due to the use of different scales 
of stellar parameters such as effective temperatures, grav- 
ities and metallicities, to the adoption of different atomic 
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parameters or the neglecting of NLTE effects. The relation 
given by the theoretical models (dotted line) is represen- 
tative of the galactic disk, since the calculated variations 
of the [O/Fe] x [Fe/H] relation from galactocentric dis- 
tances R = 4kpc to R = 14kpc are expected to be small 
(Matteucci [1996 ) . The gradient derived line (dashed line) 
was based on O/H and [Fe/H] gradients for young ob- 
jects (HII regions, hot stars and open clusters) located 
at galactocentric radii roughly in the range R = 4 kpc to 
R = 16 kpc. On the other hand, the stellar data include 
objects in the thin disk, thick disk and even some metal- 
poor halo stars. Therefore, it can be concluded that Eq. ^| 
represents fairly well the average [Fe/H] x [O/H] relation 
in the disk, so that any galactocentric variation of this 
relation is absorbed by the expected scatter, as shown in 
Fig. | 

4. The [Fe/H] x r relation 

The main difficulty in the estimate of the time variation 
of the abundance gradients from planetary nebulae lies 
in the determination of reliable ages for the central stars. 
One possibility is to use an average age-metallicity rela- 
tion, since chemical abundances from PN are relatively 
well determined, and several age-metallicity relations de- 
rived recently are similar, albeit with a considerable scat- 
ter (see Rocha-Pinto et al. 2000| for a recent discussion). 
In fact, there has been some discussion on the existence 



of such a relationship (Feltzing et al. 2001), but a critical 
analysis of the available data suggests some increase in the 
average metallicities with time, which is expected on the 
basis of the current ideas on the chemical evolution of the 
Galaxy. 

Age-metallicity relationshi ps ha ve been determin ed by 
a number of pe ople (Twarog 198C , Carlberg et al. 1985 , 
Meusinge r et al . 1991 , Edvardssone t al. 1993 Rocha-Pinto 
& Maciel |l998j Rocha-Pinto et al. |2000|), based on a va- 



riety of samples and techniques. In general, these rela- 
tionships are similar, which is particu larly true for the 
relations derived by Edvardsson et al. ( 1993 ) and the re- 
cent re sults based on chromospheric ages by Rocha-Pinto 
et al. (200C). This can be seen from Fig. 14 of Rocha- 
Pinto et al. (200C), where the new rel ation is compared 
with the relation by Edvardsson et al. ( |1993| ) adopting 1 
Gyr average bins, as well as with the results by Twarog 
( |1980D , Carlberg et al. (|1985| ) and Meusinger et al. ( |l99l| ). 
These results apply strictly to the solar neighbourhood, 
since most of the stars included in the samples are nearby 
objects with HIPPARCOS parallaxes. However, the ob- 
served scatter in these relationships is considerably large, 
amounting up to 0.26 dex for the results of Edvardsson 



et al. (1993) and about 0.13 dex for the relation obtained 
by Rocha-Pinto et al. (200C), so that it probably includes 
any differences in the corresponding relationships at dif- 
ferent galactocentric distances. This can be seen, for ex- 



ample, in Fig. 14a of Edvardsson et al. (1993), where the 
metallicity [Fe/H] is plotted against age using different 
symbols for stars at galactocentric distances R < 7 kpc, 



R= 6-7 kpc 



R= 10-11 kpc 



0.6 0.8 
log t (Gyr) 



Fig. 3. A comparison of age-metallicity relation as de- 



rived by Edvardsson et al. (1993) for stars in the galac- 
tocentric ranges 6 < i?(kpc) < 7 (top histogram) and 
10 < i?(kpc) < 11 (bottom histogram) with the estimated 
relation calculated by Eq. || (dashed lines), where we have 
used R = 6.5 kpc and R = 10.5 kpc, respectively. 



7 < i?(kpc) < 9 and R > 9 kpc, spanning radii from 4 kp c 
to 11 kpc (see also Table 14 of Edvardsson et al 



1993) 



Even though the sample at 7 < i?(kpc) < 9 is the most 
complete statistically, it can be seen that all objects can 
be reasonably fit in the average age-metallicity relation, 
in view of the large scatter. 

In order to take into account any difference in the age- 
metallicity relationship introduced by differente galacto- 
centric distances, Edvardsson et al. (1993) went a step 
further and derived an average age-metallicity-radius re- 
lation given by 



logr = 0.872 - 0.303 [Fc/H] - 0.038 R , 



(3) 



where r is the stellar age in Gyr and R is the ga lactoc entric 
distance in kpc (see eq. 19 of Edvardsson et al. 1993 ), with 
a typical uncertainty of cr(logr) ~ 0.10. This equation 
was based on stars with logr < 1.0, a condition fulfilled 
by all objects in our sample, the majority of which have 
r ~ 5 Gyr, or logr ~ 0.70. 

We can compare the ages derived by Eq. || with the 
actual average age-metallicity relation by Edvardsson et 
al. ( 1993 ) for different galactocentric distances as given 
in Table 14 of that paper. We performed this calcula- 
tion for the galactocentric radii 4 < i?(kpc) < 11, and 
Fig. H shows representative examples for 6 < i?(kpc) < 7 
and 10 < i?(kpc) < 11. The histog ram data are from 
Table 14 of Edvardsson et al. ( 1993 ) while the dashed 
lines are results from the application of Eq. ^, adopting 
R = 6.5 kpc and R — 10.5 kpc, respectively. It can be seen 
that the dashed lines produce similar results as the age- 
metallicity relation, considering that our objects have ages 
logr < 1 Gyr, peaking at r ~ 0.7 Gyr, as already men- 
tioned. Therefore, the use of Eq. || introduces a correction 
to the age-metallicity relation for the solar neighbourhood 
adjusting it to other galactocentric distances in the range 
of 4 to 11 kpc approximately. Of course, the uncertainties 
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Table 1. Results of the linear fits 



Case A 


i 
i 


1 1 

ii 


1 1 1 

in 


a log (U/ tt) 
dR 


-0.0648 


-0.0720 


-0.1164 




0.0168 


0.0057 


0.0127 


r 


-0.77 


-0.67 


-0.88 




12 


195 


27 


Case B 


I 


II 


III 


dlog(0/H) 

dR 


-0.0470 


-0.0895 


-0.0935 


°9 


0.0070 


0.0034 


0.0102 


r 


-0.64 


-0.94 


-0.75 


n 


66 


99 


69 



From the figures and the uncertainties in the slopes 
given in Table 1, it can be seen that there is a clear 
tendency for the O/H gradient to flatten out with time. 
This tendency is clear in both Cases A and B, and 
is particularly strong between Groups III and II (Case 
A) and Groups II and I (Case B). This is confirmed 
by the inspection of the uncertainties in the slopes as 
given in Table 1. For Case A, the gradients have flat- 
tened from — O.lldex/kpc to — 0.06 dex/kpc, while for 
Case B, we have —0.09 dex/kpc for the oldest group 
and —0.05 dex/kpc for the youngest one. Overall, one 
could conclude that the gradients flattened out from 
—O.lldex/kpc to —0.06 dex/kpc in about 9 Gyr, or from 
—0.08 dex/kpc to —0.06 dex/kpc in the last 5 Gyr only. 



in the derived ages are still considerably large, irrespective 
of the galactocentric radius, as can be seen by the scatter 
in the age-metallicity relations quoted here. However, in 
this study we are mainly interested in the time variation 
of the gradients, so that, in fact, only relative ages will be 
important in our analysis. Finally, as a further check of 
our results, in section 5.2. we will discuss a totally inde- 
pendent way of estimating the ages of the PN central stars 
and its effect on the time variation of the O /H abundance 
gradient. 

5. Results and discussion 

5.1. Time variation of the O/H gradient 

The PN sample has been divided into three groups of in- 
creasing ages, namely Groups I, II, and III. Since all ob- 
jects have ages under 10 Gyr, we have considered initially 
the following groups, which we will call Case A: Group I: 
< r(Gyr) < 3, Group II: 3 < r(Gyr) < 6, and Group 
III: r > 6 Gyr. 

The derived gradients are shown in Fig. |], where we 
have taken Rq — 8kpc for the galactocentric distance 
of the LSR, in a greem ent with the value adopted by 
Edvardsson et al. ( 1993 ). In fact, the main results of this 
paper are unchanged if we adopt different values such as 
the IAU recommended value Rq = 8.5 kpc, or some re- 
cently determined result such as i? = 7.6 kpc (see for ex- 
ample Maciel & Quireza 1999 ). The figure shows the sepa- 
ration of Group I (squares), II (dots) and III (crosses), and 
the linear gradients of each group (dex/kpc). The values 
of the gradients (dex/kpc), the associated uncertainties 
Cg, the correlation coefficients r and the number of ob- 
jects in each class are given in Table 1. We can see that 
Group I is underpopulated compared to the remaining 
groups, reflecting the lack of the most massive and pre- 
sumably younger objects, which affects the derived slope. 
Therefore, we have alternatively defined Case B as follows: 
Group I: < r(Gyr) < 4; Group II: 4 < r(Gyr) < 5, and 
Group III: r > 5 Gyr, so that we have approximately the 
same number of nebulae in each age group. The corre- 
sponding results are shown in Fig. || and also in Table 1. 



5.2. Ages of the central stars from their main sequence 
masses 

In view of the uncertainties involved in the estimate of 
the ages of the stars in our sample, which are basically 
due to the the adopted age-metallicity relationship, it is 
interesting to evaluate the time variation of the abun- 
dance gradients based on independent age estimates. This 
can be achieved on the basis of a correlation between the 
N/O abundances and the central star masses recently dis- 
cussed by Cazetta & Maciel ( 2000| ), which is supported 
by re cent theoretical calculations (see for example Marigo 
|2000 ). According to this relation, if N/O is the nitrogen 
abundance relative to oxygen by number of atoms, the 
core mass M c in solar masses is given by the relations 

M c ~ 0.7242 + 0.1742 log(N/0) (4) 

for log(N/0) < -0.26 and 

M c ~ 0.825 + 0.936 log(N/0) + 1.439 [log(N/0)] 2 (5) 

for log(N/0) > —0.26. These relationships are valid in the 
approximate range —1 < log(N/0) < 0.2. The main se- 
quence mass M ms can be obtained through a initial mass- 
final mass relation given by 



M C ~0.5 + 0.1M„ 



(6) 



This set of equations characterizes the so-called low-mass 
calibration (see a detailed discussion in Cazetta & Maciel 
[2000 and Maciel [2001 ) , and shows a better agreement with 
the recent discussions on the PN central star masses by 
Stasihska et al. Q1997D , Zhang ( |1993| ), with t he mass-N/O 
abundance relation by Marigo et al. (1996) and the ini- 
tial mass-final mass relations by Groenewegen & de Jong 
( |1993| ) and Blocker & Schonberner ( |1990| ). 

From Eqs. ^-^| above the main sequence mass can be 
estimated from the N/O abundances. The stellar age can 
then be derived using the average lifetimes based on stel- 
lar evolutionary models for Population I stars given by 



Bahcall & Piran (1983), that is 



\ogt(M ms ) = 10.0 - 3.6 logM ms + (logM ms ) 5 



(7) 
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Fig. 4. Radial O/H gradients from PN of Group I (squares), Group II (dots) and Group III (crosses) for Case A. 
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Fig. 5. The same as Fig. Hfor Case B. 



where t(M ms ) is given in years and M ms again in so- 
lar masses. This relation has an estimated unce rtainty of 
about 10% according to Bahcall fc Pi ran (|1983| ) or some- 



what larger according to Scalo ( 1986 ), but in any case it 
is expected to produce reliable relative ages. 

We have applied the procedure described above to the 
PN central stars in our sample, using revised N/O abun- 
dances listed by Cazetta & Maciel (200C) and Maciel & 
Chiappini (1994). Again dividing the objects into Groups 



I, II and III, we obtained the results shown in Fig. 3 
where the symbols have the same meaning as in Fig. i 
The slopes are given in the figure, and the corresponding 
uncertainties and correlation coefficients are in the range 
(T g ~ 0.013 to 0.021 and \r\ ~ 0.67 to 0.88. Although the 
available sample is smaller than in the case of the pre- 
vious figures, and some superposition can be observed in 
the different classes, it is clear that the same pattern is 
observed here, that is, the youngest objects display flatter 
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gradients, as can be seen from the slopes given on the top 
of the figure. Since the adopted procedures to estimate the 
ages in the case of Fig. [J and Fig. || are totally indepen- 
dent from each other, we have a further evidence that the 
relative behaviour of the gradients is indeed reflected by 
the results shown in Figs. ^ and g. 

5.3. Discussion 

A schematic plot of the results of Section 5.1. is shown in 
Fig. ^, where the solid lines refer to Case A and B, and the 
dotted line is a reference line at —0.07 dex/kpc. This refer- 
ence line is usually considered as an average value for the 
younger population, comprising HII regions and B stars, 
but this should be regarded with caution in view of the re- 
cent discussions on these objects as well as on open cluster 
stars and cepheids (see references in the Introduction). A 
realistic uncertainty for the slope of the young population 
can be roughly taken as 0.02 dex/kpc. For comparison 
purposes, we also include in the figur e the results of the 
theoretical models by Hou et al. ( 2000| ) based on an inside- 
out scenario for the formation of the disk with mctallicity 
dependent yields (dashed line). The agreement is remark- 
able, especially during the last 5 Gyr, which include most 
of our objects, and for which the derived ages are rel- 
atively more accurate. From these data, we can estimate 
an average flattening rate of 0.002 dex kpc -1 Gyr -1 (Case 
A) or a rate of about 0.004dex kpc -1 Gyr -1 , considering 
both cases in the last 5 Gyr, suggesting that the O/H gra- 
dient has not changed more than 30% in average during 
the last few Gyr. At earlier times, however, our results are 
consistent with a steeper rate, although the corresponding 
uncertainties are larger. 

It is interesting to compare these results with earlier 
estima tes based on planetary nebulae. Maciel & Koppcn 
( 1994 ) found some evidences for a steepening of the gra- 
dients, especially those of Ne, Ar and S, while for the 
O/H ratio some steepening could be obtained for PN of 
Peimbert types II and III only. However, the average ages 
attributed to the different PN types were very approxi- 
mate, and no effort was made to establish individual ages. 
These av erage values have also been used by Maciel & 
Quircza ( 1999 ), leading to the conclusion that the gra- 
dients steepen out in time at a ratio of about 0.004 dex 
kpc -1 Gyr" 1 in the Galaxy, a result which is not con- 
firmed by the present paper. The main difference between 
the present work and these early attempts, apart from our 
larger and more accurate sample, refers to the definition 
of the age groups. Maciel and Koppen (1994) adopted lit- 
erally the Peimbert classification scheme, which, strictly 
speaking, attributes different and increasing ages for the 
progenitor stars of PN of types I, II and III, respectively. 
In this work, we have not assumed any correlation between 
the PN type and the progenitor age, and made an effort 
to derive individual ages, so that we expect the present re- 
sults to be more reliable. We notice that most Type I PN 
in our sample belong to age Group II, most Type II PN 



are spread between Groups I and II, while PN of Type III 
belong either to Groups II or III. Finally, most of the ob- 
jects in Group III are indeed Type III PN (Case A). We 
can then conclude that within the disk PN of a given type 
there are objects with a reasonable age span. In other 
words, there seems to be some overlapping in the progeni- 
tor masses - and ages - within the Peimbert types, which 
was not accounted for by the earlier attempts to estimate 
the time variation of the gradients, and probably explains 
the different results presented here as compared with the 
earlier work. In fact, as recently discussed by Peimbert & 



Carigi (1998), some overlapping occurs in the masses of 
the progenitors of PN of types II and III, which supports 
our present views. 
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ages of the central stars have been calculated as described in Section 5.2. 
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